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What do WMAP and SDSS really tell about inflation? 
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We derive new constraints on the Hubble function H(cf>) and subsequently on the inflationary 
potential V(4>) from WMAP 3-year data combined with the Sloan Luminous Red Galaxy sur- 
vey (SDSS-LRG), using a new methodology which appears to be more generic, conservative and 
model-independent than in most of the recent literature, since it depends neither on the slow-roll 
approximation for computing the primordial spectra, nor on any extrapolation scheme for the po- 
tential beyond the observable e-fold range, nor on additional assumptions about initial conditions 
for the infiaton velocity. This last feature represents the main improvement of this work, and is 
made possible by the reconstruction of H((f>) prior to V((j>). Our results only rely on the assumption 
that within the observable range, corresponding to ~ 10 e-folds, inflation is not interrupted and the 
function H((j>) is smooth enough for being Taylor-expanded at order one, two or three. We conclude 
that the variety of potentials allowed by the data is still large. However, it is clear that the first 
two slow-roll parameters are really small while the validity of the slow-roll expansion beyond them 
is not established. 

PACS numbers: 98.80.Cq 
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Cosmic inflation was introduced as a simple and aes- 
thetically elegant scenario of the early Universe evolu- 
tion which is capable of explaining its main properties 
observed at the present time 0, 0, y|, 0, H, @] ■ As a very 
important byproduct it provides a successful mechanism 
for the quantum-gravitational generation of primordial 
scalar (density) perturbations and gravitational waves 



Refs. [27|, |2a |29| for specific inflationary models). The 
extrapolation over the full duration of inflation is more 
constraining than the data alone. Instead, Ref. [13] fo- 
cused only on the observable part of the potential to see 
up to what extent current data really constrains inflation. 

For this class of models, the evolution of a spatially flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) uni- 



@, H, Hi, E3, Hi 12, EH- The Fourier power spectrum verse can be described by [3(3, 31 1 



Viz{k) of the former ones is observed today in the cos- 
mic microwave background (CMB) and the large scale 
structure (LSS). Vice versa, at present the CMB and 
the LSS provide the only quantifiable observables which 
can confirm or falsify inflationary predictions. That is 
why matching concrete inflationary models to observa- 
tions has become one of the leading quests in cosmology. 

In the simplest class of inflationary models, inflation is 
driven by a single scalar field cf> (an infiaton) with some 
potential V(4>) which is minimally coupled to the Ein- 
stein gravity. For these models, some new conservative 
bounds on V{<p) were presented recently in [TiJ. Until 
then, most post-WMAP3 studies concerning V(4>) relied 
on the slow-roll approximation in the calculation of per- 
turbation power spectra and their relation to values of 
<f> during inflation [h| [~ 
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or made an extrapolation of V(6) from the observable 
window till the end of inflation 24, 25|, 26| (a numerical 
integration of exact wave equations for perturbations to 
obtain primordial power spectra was also permormed in 
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whenever </> ^ and not specifically during inflation (so 
H'(<f>) ^ 0, too). Here H(<p(t)) = a/a, a(t) is the FLRW 
scale factor, a dot denotes the derivative with respect to 
the cosmic time t, a prime with respect to an argument, 
and we have set Gm 2 P = h = c = 1. If V(<f>) is consid- 
ered as the defining quantity, the initial conditions for 
generating the observable window are determined by the 
set {(/>ini, V{4>)}. In Ref. the infiaton potential was 
parametrized as a Taylor expansion up to some order, to 
see up to what extent the potential can be constrained 
by pure observations. However, in order to reduce the 
number of free parameters, 0j n i was fixed for each model 
by demanding that the infiaton follows its attractor so- 
lution just when the observable modes exit the horizon. 
In practice this means that the results of Ref. [3] as- 
sumed that inflation started at least a few e-folds before 
the observable modes left the horizon. These precurring 
e-folds led to a slightly stronger bound on the potentials 
than the data itself could actually give, although this 
extra constraining power stands in no proportion to an 
extrapolation over the full duration of inflation. 
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Eqs. (JT][2j) however show that when one considers H(4>) 
as the defining quantity, all initial conditions arc already 
uniquely set by H (</>). Moreover, the slow-roll conditions 
which require, in particular, that the first term in the 
rhs of Eq. ^ is much less than the last one need not be 
imposed ab initio. In this Letter we derive the bounds on 
H((f>) during observable inflation using its Taylor expan- 
sion at various orders. We infer for this some constraints 
on V((f> ) un der an even more conservative approach than 
in Ref. [lj| , since the present method requires absolutely 
no extrapolation outside of the observable region (cither 
forward or backward in time). Our only restriction is to 
assume that observable cosmological perturbations orig- 
inate from the quantum fluctuations of a single inflaton 
field, which dynamics during observable inflation is com- 
patible with a smooth, featureless H{4>). 

Method. We used the publicly available code COS- 
momc [H to do a Monte Carlo Markov Chain (MCMC) 
simulation. We added a new module (released at 
|http : / /wwwlapp . in2p3 . f r/~va lkenbu/inf lationH/) 
which computes numerically the primordial spectrum of 
scalar and tensor perturbations for each given function 
H((j) — 0»), where 0* is an arbitrary pivot scale in field 
space. This module is simpler than the one in Ref. [l4j |. 
since the code never needs to find an attractor solution 
of the form 0(0). The comoving pivot wavenumber 
is fixed once and for all to be fc* = 0.01 Mpc -1 , 
roughly in the middle of the observable range. Pri- 
mordial power spectra are computed in the range 
[fcmin,fcmax] = [5 x 10~ 6 , 5] Mpc -1 needed by CAMB, 
imposing that fc* leaves the Hubble radius when = 0*. 
In practice, this just means that for each model the code 
normalizes the scale factor to the value a* = fc*/-ff* 
when — 0*. Note that by mapping a window of infla- 
tion to a window of observations today, our approach 
is independent of the mechanism of reheating. The 
evolution of each scalar/tensor mode is given by 
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with r\ = J dt/a(t) and z$ = acj)/H for scalars, zt = a 
for tensors. The code integrates this equation start- 
ing from the initial condition £s.t = e _?fe, '/v / 2fc when 
k/aH = 50, and stops when the expression for the ob- 
served scalar/tensor power spectrum freezes out in the 
long-wavelength regime. More precisely, the spectra are 
given by 
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and integration stops when [dlnVn^/dlna] < 10~ 3 . If 
for a given function H{(f> — 0*) the product aH can- 
not grow enough for fullfilling the above conditions, the 
model is rejected. In addition, we impose that aH grows 
monotonically, which is equivalent to saying that inflation 
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TABLE I: Bayesian 68% confidence limits for ACDM infla- 
tionary models with a Taylor expansion of H(<f) — 0„) at order 
n — 1,2,3 (with the primordial spectra computed numeri- 
cally). The last line shows the maximum likelihood. The first 
four parameters have standard definitions (see e.g. [3|)- 



is not interrupted during the observable range. If these 
conditions are satisfied, the power spectra are compared 
to observations. 

We choose to parametrize H as a Taylor expansion 
with respect to 0—0* up to a given order n vary- 
ing between one and three (this choice of background 
parametrization is equivalent to that in Ref. [ls| , as long 
as no extrapolation is made). Note that for n > 1 
such an assumption excludes and H' becoming zero 
at some value = 0i in the range involved since then 
i?(0) would acquire a non-analytic part beginning from 
the term oc |0 — 0i | 3 / 2 (with L(0) being totally analytic 
at this point) [4l|. As a cosmological background we 
used the standard ACDM-model with the free parame- 
ters shown in Table U 

Results for H(<fi — 0*). In Fig. Q] we show the prob- 
ability distribution of each parameter marginalized over 
the other parameters. The corresponding 68% confidence 
limits are displayed in Table [H as well as the minimum 
of the effective x 2 f° r each model. This minimum does 
not decrease significantly when n increases, which reflects 
the fact that current data are compatible with the sim- 
plest spectra and potentials, but derivatives up to H'" 
can be constrained with good accuracy. Note that it 
would be very difficult to give bounds directly on the 
set {H,H',H",H"',...}: indeed, these parameters are 
strongly correlated by the data, because physical effects 
in the power spectra depend on combinations of them. 
For example, at the pivot scale, the scalar amplitude 
is mainly determined by (H 2 /Hl) 2 and the tensor-to- 
scalar ratio r = Vh/Vn by (if* /if*) 2 . The scalar tilt 
ns further depends on H'J/H«, and the scalar running 
on H"' H^/ H 2 . The Markov Chains can converge in a 
reasonable amount of time only if the basis of parame- 
ters (receiving flat priors) consists in functions of each 
of the above quantities, or linear combinations of them. 
However, we also show in the last plot of Fig.[T]thc distri- 
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FIG. 1: Probability distribution for the eight independent 
parameters of the models considered here, normalized to a 
common arbitrary value of P max . The ninth plot shows a re- 
lated parameter (with non-flat prior) : namely, the value of the 
expansion rate when the pivot scale leaves the horizon during 
inflation. Our three runs n = 1,2,3 correspond respectively 
to the dashed red, dotted blue and solid black lines. The 
data consists of the WMAP 3-year results QJ, H, H, HI and 
the SDSS LRG spectrum [3(| . The first four parameters have 
standard definitions (see e.g. [3), and \nA is a shortcut 
notation for the parameter defined in the fifth line of Table I. 



bution of H*: this information is useful since the energy 
scale of inflation is given by A = (SH^n^p/Sn) 1 / 4 , but 
the displayed probability should be interpreted with care 
since this parameter has a non-flat prior. 

The run n = 1 is not very interesting. Indeed, impos- 
ing H" and higher derivatives to vanish leads to a one- 
to-one correspondence (at least in the slow-roll limit) be- 
tween the amplitude and the tilt of the scalar spectrum. 
This feature is rather artificial and unmotivated. It ex- 
plains anyway why the parameter has exceptionally 
a lower bound in the n = 1 case 42J. Much more inter- 
esting is the n = 2 case for which the tensor ratio, scalar 
amplitude and scalar tilt are completely independent of 
each other, and the n = 3 case for which even the tilt 
running has complete freedom. The runs for n = 2 and 
?i = 3 nicely converged and constitute the main result 
of this work. Note also that the middle-right and lower- 
right graphs in Fig. 1 are compatible with each other 
in the following sense: though H'* may not reach zero 
under our assumption, the quantity H' * /H* may be 



arbitrarily small if H* is allowed to be arbitrarily small, 
too. Thus, for cases n= 2,3 when H* is not suppressed 
at zero argument, H' * /H* is not suppressed there, too. 

The probability distribution for combinations of H*, 
H* and H" are robust in the sense that they do not 
change significantly when one extra free parameter H"' 
is included: this indicates that they are directly con- 
strained by the data. We tried to include an additional 
parameter (H'J" / H*)(H^/ 'H*) 2 m P , but then our Markov 
Chains did not converge even after accumulating of the 
order of 10 5 samples. We conclude that current data do 
not have the sensitivity required to constrain H(4>) be- 
yond its third derivative and to establish the validity of 
the slow-roll approximation beginning from this order. 
On the other hand, the first two slow-roll parameters 
e(4>) = H' 2 m%/4:TrH 2 and rj(4>) = H"mpj^H are really 
small over the observed range (tilde is used here to avoid 
mixing with the conformal time 77). The next parame- 
ter £ = 2 \ H = H'"H'mp/(4ir) 2 H 2 is also small, ~ 0.01, 
though being of the order of e and \f]\, not e 2 or fj 2 as 
would follow from the standard slow-roll expansion. This 
smallness explains why our results for these parameters 
are similar to those obtained for the same background 
H((f>) but using the slow-roll approximation to calculate 
the power spectra [2l| (and to those in [24[ , too) although 
some important differences exist. 

Results for V{4> — 0*). We further processed our 
n = 1,2,3 runs in order to reconstruct the inflaton 
potential. For each run, we kept only 68% or 95% of 
the models with the best likelihood, and computed the 
corresponding inflaton potentials using Eq. ([2]). Note 
that the problem is fully symmetric under the reflec- 
tion (0—0*) <-> —((/> — (/>»). We choose to focus on one 
half of the solutions, corresponding to <fi > and hence 
VI > 0. Our results are shown in Fig. [H They appear 
to be compatible with those of Ref . [14[ , although a de- 
tailed comparison is difficult: first, the current method 
is more conservative, and second, a given order in the 
Taylor-expansion of H(<fi — </>*) is not equivalent to an- 
other order in that of V{<j> — 0*). Our results are also 
difficult to compare with those of Ref. (26[, since these 
authors choose to present their full allowed potentials ex- 
trapolated till the end of inflation: in principle, our Fig. [2] 
can be seen as a zoom on the directly constrained, small 
<t> region in their Fig. 2. 

Our results could give the wrong impression that all 
preferred potentials are concave. This comes from the 
fact that in the representation of Fig. [5J many interest- 
ing potentials are hidden, since they almost reduce to the 
point (K,A0) — ► (0,0). Indeed, as long as the tensor- 
to-scalar ratio is not bounded from below, many low- 
energy inflationary models with very small and HI 
(and hence tiny variation of the inflaton field during the 
observable e-folds) are perfectly compatible with obser- 
vations. It is straightforward to show that models leading 
to ns < 1 and small r correspond to convex potentials 
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FIG. 2: Allowed inflationary potentials V((j> — <f>*) inferred 
from each of our n = 1,2,3 runs. For each case, the light 
colour corresponds to models allowed at the 68% confidence 
levels, and the dark colour to the 95% level. The inner (bluish) 
region is obtained for n = 1, the intermediate (greenish) one 
for n — 2 and the outer (reddish) one for n — 3. Each po- 
tential is plotted between the two values <j>i and <p2 corre- 
sponding to Hubble exit for the limits of the observable range 
[fei,fc2] = [2 x 10 _4 ,0.1] Mpc -1 : so we only see here the ac- 
tual observable part of each potential. Note that this figure 
shows only one half of the possible solutions: the other half 
is obtained by reflection around <j> = <f>* . 



(like e.g. new inflation with V = Vq — \<fi n , or one- loop 
hybrid inflation with V = Vq + A In <fi) , while models with 
same ns and larger r derive from concave potentials (like 
e.g. monomial inflation V = \(f) a ). Current data favor 
ns < 1, and the upper bound on r is too loose for differ- 
entiating between these two situations. So, our allowed 
potentials can be split in two subsets: low-energy convex 
potentials and high-energy concave potentials, as illus- 
trated in Fig. [3l in which we rescaled all allowed poten- 
tials to the same variation in V and <j>. More generally, 
this large degeneracy in potential reconstruction reflects 
the fact that an infinitely precise measurement of the 
scalar spectrum V-ji would only constrain the function 
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(in the slow-roll approximation). This is not sufficient 
for inferring the correspondence between k and <fr, and 
hence for a unique determination of H(<fr) and V(<f>). It 
is necessary to measure also the tensor spectrum, equal 
to 
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in the same approximation, in order to diminish this de- 
generacy (see the related discussion in Ref. [3lJ )• In the 
slow- roll approximation, the knowledge of Vhik) leads 



FIG. 3: Allowed inflationary potentials V(4> — </>*) with the 
same colour/shade code as in Fig. [2] but a different choice of 
axes: each potential is now rescaled to the same variation in 
V an <f> space. This shows that many allowed potentials are 
actually convex. The outer region still corresponds n = 3, the 
intermediate one to n — 2 and the inner (quasi-linear) one to 
n = 1. 



to the unambiguous determination of H((f>). However, 
the question how unique the determination of H{<p) is, 
even from both Vu{k) and Vh{k) in the generic case be- 
yond slow-roll, is still open because of the existence of 
many -/?(</>) leading to the same perturbation spectra 
which may not be obtained from the slow-roll expansion 
at all |38| ■ Still, since the difference of these additional 
solutions from slow-roll ones is, in some sense, exponen- 
tially small for small slow-roll parameters, their existence 
might appear not significant from the observational point 
of view. 
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